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Abstract 

Yeast glycolysis is considered the prototype of dissipative biochemical oscillators. In cellular 
conditions, under sinusoidal source of glucose, the activity of glycolytic enzymes can display either 
periodic, quasiperiodic or chaotic behavior. In order to quantify the functional connectivity for the 
glycolytic enzymes in dissipative conditions we have analyzed different catalytic patterns using the 
non-linear statistical tool of Transfer Entropy. The data were obtained by means of a yeast glycolytic 
model formed by three delay differential equations where the enzymatic speed functions of the irre- 
versible stages have been explicitly considered. These enzymatic activity functions were previously 
modeled and tested experimentally by other different groups. In agreement with experimental con- 
ditions, the studied time series corresponded to a quasi-periodic route to chaos. The results of the 
analysis are three-fold: first, in addition to the classical topological structure characterized by the 
specific location of enzymes, substrates, products and feedback regulatory metabolites, an effective 
functional structure emerges in the modeled glycolytic system, which is dynamical and character- 
ized by notable variations of the functional interactions. Second, the dynamical structure exhibits a 
metabolic invariant which constrains the functional attributes of the enzymes. Finally, in accordance 
with the classical biochemical studies, our numerical analysis reveals in a quantitative manner that 
the enzyme phosphofructokinase is the key-core of the metabolic system, behaving for all conditions 
as the main source of the effective causal flows in yeast glycolysis. 
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Author Summary 

The understanding of the effective functionality that governs the enzymatic self-organized processes in 
cellular conditions is a crucial topic in the post-gcnomic era. A number of measures have been proposed 
for the functionality and correlations between biochemical time series. However, functional correlations 
do not imply effective connectivity and most synchronization measures do not distinguish between causal 
and non-causal interactions. In recent studies, Transfer Entropy (TE) has been proposed as a rigorous, 
robust and self-consistent method for the causal quantification of the functional information flow among 
nonlinear processes. Here, we have used TE to establish the effective functional connectivity of yeast 
glycolysis under dissipative conditions. Concretely, we have applied this method for a quantification of 
how much the temporal evolution of the activity of one enzyme helps to improve the future prediction of 
another. In the enzymatic activities, the oscillatory patterns of the metabolic products might have causal 
information which can be appropriately read-out by the TE. We have performed numerical studies of yeast 
glycolysis under dissipative conditions and found the emergence of a new kind of dynamical functional 
structure, characterized by changing connectivity flows and a metabolic invariant that constrains the 
activity of the irreversible enzymes. 
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Introduction 

Yeast glycolysis is one of the most studied dissipative pathways of the cell; it was the first metabolic 
system in which spontaneous oscillations were observed [IJ[2] , and the study of these rhythms allowed the 
construction of the first dynamic model where the kinetics of an enzyme was explicitly considered [2l|4]. 
More concretely, the main instability-generating mechanism in the yeast glycolysis is based on the self- 
catalytic regulation of the enzyme phosphofructokinase (31E3I6] . 

Glycolysis is the central pathway of glucose degradation which is implied in relevant metabolic pro- 
cesses, such as the maintenance of cellular redox states, the provision of ATP for membrane pumps and 
protein phosphorylation, biosynthesis, etc; and its activity is linked to a high variety of important cellular 
processes, e.g., glycolysis has a long history in cancer cell biology [7] and cell proliferation [5], there is a 
correlation between brain aerobic glycolysis and Amyloid-/? plaque deposition which might precede the 
clinical manifestations of the Alzheimer disease [5] , the glycolytic inhibition abrogates epileptogenesis [TU] , 
and glycolysis is also related with oxidative stress QT] and apoptosis |12| . 

Over the last 30 years a large number of different studies focused on different molecular mechanisms 
allowing for the emergency of self-organized glycolytic patterns [T5HT5] . Nevertheless, despite the intense 
advance of the knowledge of these metabolic structure, we still lack a quantitative description in cellular 
conditions of the effective functional structure and the causal effects among the enzymes. 

In this paper, to go a next step further in the understanding of the relationship between the classical 
topological structure and functionality we have analyzed the effective connectivity of yeast glycolysis, 
which in inter-enzyme interactions accounts for the influence that the activity of one enzyme has on the 
future of another [T§1 - |2"3"] . 

For this purpose, we considered a yeast glycolytic model described by a system of three delay- 
differential equations in which there is an explicit consideration of the speed functions of the three 
irreversible enzymes hexokinase, phosphofructokinase and pyruvatekinase. These enzymatic activity 
functions were previously modeled and tested experimentally by other different groups }3"l l24l[2"5] . 

We have obtained time scries of enzymatic activity under different sources of the glucose input flux. 
The data corresponded to a typical quasi-periodic route to chaos which is in agreement with experi- 
mental conditions |26j . The dynamics of the glycolytic system changes substantially trough this route, 
which allows for a better comparison of the enzymatic processes in periodic, quasi-periodic and chaotic 
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conditions. 

Using the non- linear analysis techniques such as Transfer Entropy and Mutual Information |28j . 
we have analyzed the glycolytic series and quantified the effective connectivity of the enzymes. 

The results show that in the numerical analysis of yeast glycolysis, under dissipative conditions, 
a effective functional structure emerges which is characterized by changing connectivity flows and a 
metabolic invariant that constraints the activity of the irreversible enzymes. 

Results 

The monitoring of the fluorescence of NADH in glycolyzing baker's yeast under sinusoidal glucose input 
flux, have shown that quasi-periodic time patterns are common at low amplitudes of the input and for 
high amplitudes chaotic behaviours emerge f2"W3"0"] . 

In order to simulate these metabolic processes, the system is considered under periodic input flux with 
a sinusoidal source of glucose S = So + Asin(wt). Assuming the experimental value of So = 6mM/h [31 j . 
after dividing by K m 2 (the Michaelis constant of phosphofructokinase, see for more details Materials and 
Methods) we have obtained the normalized input flux So = 0.033 Hz. 

Under these conditions, a wide range of different types of dynamic patterns can emerge as a function 
of the control parameter, hereafter the amplitude A of the sinusoidal glucose input flux ■ In 

particular, it is observed a quasi-periodic route to chaos (cf. left panel in Fig. [5]); thus for A= 0.001 the 
biochemical oscillator exhibits a periodic pattern (Figure 2a). An increment of the amplitude to A= 0.005 
provokes a Hopf bifurcation generating another fundamental frequency, as a consequence, quasi-periodic 
behaviors emerge (Figure 2b). Above A= 0.021, complex quasi-periodic oscillations appear (Figure 
2c). After a new Hopf bifurcation the originated dynamical behavior is not particularly stable and small 
perturbations produce deterministic chaos (A= 0.023, Figure 2d), as predicted by Ruelle and Takens |34) . 
This route is in agreement with experimental conditions |26j . 

To go a next step further in the understanding of the relationship between the classical topological 
structure and effective functionality we have analyzed by means of non-linear statistical tools the catalytic 
patterns belonging to this scenario to chaos, and for each transition represented in the Figure 2 we have 
obtained three time series corresponding to the variables a, j3 and 7 (12 in total), which denote respectively 
the normalized concentrations of glucose-6-phosphate, fructose 1-6-bisphosphate and pyruvate. 
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Effective functionality 

Transfer Entropy (TE) quantifies the reduction in uncertainty that one variable has on its own future 
when adding another. This measure allows for a calculation of the functional influence in terms of effective 
connectivity between two variables [57| . The analysis of the glycolytic data by means of the TE method 
are shown in Table I. The 4D vectors in square brackets correspond to the results obtained for the 4 
different amplitudes of the considered glucose input flux, A=[0.001;0.005;0.021;0.023]. 

The values of functional influence are ranging in 0.58 < TE < 1.00, with mean= 0.79 and standard 
deviation=0.12, what indicates in general terms a high effective connectivity in the enzymatic system. 
The minimum value 0.58 corresponded to the causality flow between E3 and E2 when a simple periodic 
behavior emerges. However, the functional connectivity from E2 to E3 shows the maximum value, 
achieved in all considered conditions of the glucose input flux. 

The glycolytic effective connectivity is illustrated in the right panel of Fig. 2. The arrows width is 
proportional to the TE between pairs of enzymes. The values change trough the quasi-periodic route to 
chaos, remarked from E3 to E2 by black dashed circles, [0.58;0.84;0.61;0.66]. 

In all cases analyzed, the values of TE present a maximum statistical significance (pvalue=0). 

Total Information flows and the functional invariant 

Next, we have measured the total information flow, defined as the total outward of Transfer Entropy 
arriving to one enzyme minus the total inward. Positive values mean that that enzyme is a source of 
causality flow and negative flows are interpreted as sinks or targets. The results of the total information 
flows are shown in Table II (pvaluc=0). The maximum source of total transfer information (0.41) corre- 
sponds to the E2 enzyme (phosphofructokinase) for A=0.021, when complex quasi-periodic oscillations 
appear in the glycolytic system. 

For all conditions the enzyme E2 (phosphofructokinase) is the main source of effective influence and 
the enzyme E3 (pyruvatekinase) a sink, which could be interpreted as a target from a point of view based 
on its effective functionality. The enzyme Ei (hexokinase) is less constrained, and it has a flow close to 
zero for all conditions. 

The attributed role to each enzyme, namely E2 the source, E3 the sink and Ei no-constrained is an 
invariant and preserved trough the whole route to chaos. 



G 



Functional Synchronization 

Time correlations allows for quantification about how much two time series are statistically independent. 
According to that, we have measured the time pairwise correlations in the enzymatic system, and the 
corresponding results are shown in Table III. The main finding is that E 2 and E 3 are highly synchronized 
(correlation=0.90, pvalue=0) and Ei is anti-synchronizcd with both E2 and E3 (respectively, correlation 
equals -0.65 and -0.66, pvalue=0). 

These values of time correlations were almost constant trough the quasi-periodic route to chaos and 
established that the activities of E2 and E3 are grouped to the same function, being activated at similar 
time and oppositely to Ei. 

Redundancy and uncertainty reduction 

The Mutual Information (MI) quantifies how much the knowledge of one variable reduces the entropy or 
uncertainty of the another [55] . The analysis of the glycolytic data by means of this method are shown 
in Table IV. 

The high values of MI (close to 0.50) proved a high informative redundancy between the pairs of 
enzymes. So, the number of bits of information transferred from one enzyme to another is much larger 
than the actually needed. 

The values in the principal diagonal of Table IV represent the uncertainty for each variable. We have 
found these values gradually descending, H(E 1 )=[1.00;1.00;1.00;1.00], H(E 2 )=[0.85;0.84;0.85;0.86] and 
H(E2)= [0.76;0.74;0.76;0.78], which is indicative of the uncertainty in the enzymatic activity patterns 
belonging to Ei, E2 and E3 is reduced monotonously for all analyzed conditions. 

The values of MI have a maximum statistical significance (pvalue=0). 

Finally, we have computed the Mutual Information between the glucose input fluxes and the activity 
patterns of the different enzymes. In all cases, the MI was equal to zero, proving that the oscillations 
of the glucose were statistical independent of the glucose- 6-phosphate, fructose 1-6-biphosphate and 
pyruvate, products of the main irreversible enzymes of glycolysis. 
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Discussion 

In this paper we have quantified essential aspects of the effective functional connectivity among the main 
glycolytic enzymes in dissipativc conditions. 

First, we have computed under different source of glucose the causality flows in the metabolic system. 
This level of the functional influence accounts for the contribution of each enzyme to the generation of 
the different catalytic behavior and adds a directionality in the influence interactions between enzymes. 

The results show that the flows of functional connectivity change significantly during the different 
metabolic transitions analyzed, exhibiting high values of transfer entropy, and in all considered cases, the 
enzyme phosphofructokinase (E2) is the main source of effective causality flow; the pyruvatckinase (E3) 
is the main sink of information flow; the hexokinase (Ei) has a quasi-zero flow, meaning that, the total 
information arriving to Ei goes out to either E2 or E3. 

The maximum source of total transfer information (0.41) corresponds to the E2 enzyme (phospho- 
fructokinase) at the edge of chaos, when complex quasi-periodic oscillations emerge (cf. Fig. 2). This 
finding seems to be consistent with other studies which show that when a dynamical system operates in 
the frontier between order (periodic behavior) and chaos its complexity is maximal |35U36| . 

The level of influence in terms of causal interactions between the enzymes is not always the same 
but varies depending on the substrate fluxes and the dynamic characteristics emerging in the system. 
In addition to the glycolitic topological structure characterized by the specific location of enzymes, sub- 
strates, products and regulatory metabolites there is an functional structure of information flows which 
is dynamic and exhibit notable variations of the causal interactions. 

Another aspect of the glycolitic functionality was observed during the quantification of the Mutual 
Information, which measures how much the uncertainty about the one enzyme is reduced by knowing the 
other; we found that the uncertainty for Ei, E2 and E3 monotonously decreased for all the values of the 
periodic glucose input-flux. 

Second, the numerical results show that for all analyzed cases the maximum effective connectivity 
corresponds to the Transfer Entropy from E2 to E3, indicating the biggest information flow in the multi- 
enzyme instability-generating system. This is also corroborated by the measure of correlation between 
the different pairs of series which shows that E2 and E3 are highly correlated, or synchronized (corre- 
lation^. 90 pvalue=0) and Ei is anti-correlated with both E2 and E3 (respectively, correlation=-0.65 
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pvaluc=0 and corrclation=-0.66 pvaluc=0). The values of time correlations establish that the activities 
of E 2 and E 3 are grouped to the same function, being activated at similar time and oppositely to Ex- 

Third, our analysis allows for a hierarchical classification in terms of what glycolytic enzyme is im- 
proving the future prediction of what others, and the results reveals in a quantitative manner that the 
enzyme E2 (phosphofructokinase) is the major source of causal information and represents the key-core 
of glycolysis. The second in importance is the E3 (pyruvatckinase). 

From the biochemical point of view the E2 (phosphofructokinase) has been commonly considered as 
a major checkpoint in the control of glycolysis [37ll38| . The main reason for this generalized belief is that 
this enzyme exhibits a complex regulatory behavior that reflects its capacity to integrate many different 
signals |39j ; from a dissipativc point of view, this enzyme catalyzes a reaction very far from equilibrium 
and its self-catalytic regulation it has been considered the main instability-generating mechanism for 
the emergence of oscillatory patters in glycolysis [5J. The functional studies presented here confirm in a 
quantitative manner that the E2 (phosphofructokinase) is the key-core of the pathway, and our results 
make stronger and expand the classical biochemical studies of glycolysis. 

Forth, the dynamics of the glycolytic system changes substantially trough the quasi-periodic route to 
chaos when the amplitude of the input-flux varies. However, the hierarchy obtained by transfer entropy, 
E2 the flow, E3 the sink and Ei a quasi-zero flow, is preserved during this route and seems to be an 
invariant. This functional invariant of a metabolic process may be important for the understanding of 
functional enzymatic constraints in cellular conditions; but this issue requires other additional studies. 

Finally, we want to emphasize that Transfer Entropy as a quantitative measure of effective causal 
connectivity can be a very useful tool in studies of enzymatic processes that operate far from equilibrium 
conditions. Moreover, many experimental observations have shown that the oscillations in the enzymatic 
activity seem to represent one of the most striking manifestations of the metabolic dynamic behaviors, 
of not only qualitative but also quantitative importance in cells (further details in Appendix III). 

Transfer Entropy is able to detect the directed exchange of causality flows among the irreversible 
enzymes which might allow for a rigorous quantification of the effective functional connectivity of many 
dissipative metabolic processes in both normal and pathological cellular conditions. 

The TE method applied to our numerical studies of yeast glycolisis shows the emergence of a new kind 
of dynamical functional structure which is characterized by changing connectivity flows and a metabolic 
invariant that constrains the activity of the irreversible enzymes. 
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The understanding the effective connectivity of the metabolic dissipativc structures is crucial to ad- 
dress the functional dynamics of cellular life. 

Methods 
Model 

In Fig. [T]arc represented the main enzymatic processes of yeast glycolysis (the irreversible stages) with the 
enzymes arranged in series. When the metabolite S (glucose) feeds the system, it is transformed by the 
first enzyme Ei (hcxokinasc) into the product Pi (glucosc-6-phosphate). The enzymes E2 (phosphofruc- 
tokinase) and E3 (pyruvatekinase) are allosteric, and transform the substrates P'i (fructose 6-phosphate) 
and P' 2 (phosphocnolpyruvate) in the products P2 (fructose 1-6-bisphosphate) and P3 (pyruvate), respec- 
tively. The step P2 — > P' 2 represents reversible activity processes, reflected in the dynamic system by the 
functional variable /?'. A part of Pi does not continue in the metabolic system, and is removed with a 
rate constant of qi which is related with the activity of pentose phosphate pathway; likewise, q2 is the 
rate constant for the sink of the product P3 which is related with the activity of pyruvate dehydrogenase 
complex. 

The main instability-generating mechanism in yeast glycolysis is the self-catalytic regulation of the 
enzyme E2 (phosphofructokinase) , specifically, the positive feed-back exerted by the reaction products, the 
ADP and fructose- 1,6-bisphosphate [511516] . From a strictly biochemical point of view, E2 is also considered 
the main regulator enzyme of glycolysis |39| . The second irreversible stage for its regulatory importance 
is catalyzed by the enzyme E3 (pyruvatekinase) which is inhibited by the ATP reaction product |39j . 
Finally, the third irreversible process corresponds to the first stage the enzyme Ei (hcxokinasc) which is 
dependent on the ATP. 

In the determination of the enzymatic kinetics of the enzyme Ei (hexokinase) the equation of the 
reaction speed dependent on glucose and ATP has been used [23]. The speed function of the allosteric 
enzyme E2 (phosphofructokinase) was developed in the framework of the concerted transition theory [3] . 
The reaction speed of the enzyme E3 (pyruvatekinase), dependent on ATP and phospoenolpyruvate, was 
also constructed on the allosteric model of the concerted transition [25] . 

To study the kinetics of the dissipative glycolytic system we have considered normalized concentra- 
tions; a, /3 and 7 denoted respectively the normalized concentrations of Pi, P2 and P3. For a spatially 
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homogeneous system the time-evolution is described by the following three delay differential equations: 

da 

— = zio-i0i(ju) - <7 2 02(a,/3) - qia 

— = z 3 cr 3 3 (/3,/3',/i) - q 2 7 (1) 

where the functional variables /?' and /z reflect the normalized concentrations of P' 2 (phosphoenolpyruvate) 
and ATP respectively. The three main enzymatic functions are the following: 



<t>\(p) 



^SK d3 



(K 3 K 2 + /xK m iK d3 + SK 2 + ^SKda) 

q(l + q)(l + di/3) 2 
Li(l + ca) 2 + (l + a) 2 (l + di^) 2 



A, (R R> \ d 2 /3'(l + d 2 /?') 3 

</>3{P,P,M = — — ——4 71 —^4 ( 2 ) 

L 2 (l + d 3 A*) +(l + d 2 /3) 



and 



P = W-Ai)) 

M = hG8(t-A 2 )). (3) 

The constants <7i, tr 2 and cr 3 correspond to the maximum activity of Ei, E 2 and E 3 (Kni ; Vm2 and Kn3) 
divided by the Michaclis constants of each enzyme, respectively K m i, K m2 and K m3 . The constants z's 
arc defined as zi = K m i/K m2 , z 2 = K m2 /K m3 and z 3 = K m3 /Kd 3 , with Kd 3 representing the dissociation 
constant of P 2 by E 3 . The constants d's are di = K m3 /Kd 2 , d 2 = K m3 /Kd 3 and d 3 = Kd 3 /Kd4, with 
Kd4 representing the dissociation constant of ATP; Li and L 2 are respectively the allosteric constant of 
E 2 and E 3 ; c is the non-exclusive binding coefficient of the substrate Pi. More details about parameter 
values and experimental references are given in Appendix I. 

From the dissipative point of view the essential enzymatic stages are those that correspond to the 
biochemical irreversible processes |40j and to simplify the model, we did not consider the intermediate 
part of glycolysis belonging to the enzymatic reversible stages. In this way, the functions f and h are 
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supposed to be the identity function. Thus, 



fx = 7(t - A 2 ) 



(4) 



The initial functions present a simple harmonic oscillation in the following form: 



a {t) 



A + Bsm(2ir/P) 



C + D sin(27T / P) 



70 W 



E + Fsm(2ir/P) 



(5) 



with A = 26, B = 12, C = 12, L> = 10, E = 7, F = 6 and P = 534. 

The dependent variables a, /3 and 7 were normalized dividing them by K m 2, K m 3 and Kd3, and 
the parameters Ai and A2 are time delays affecting the independent variable (see for more details the 
Appendix II). 

The numerical integration of the system was performed with the package ODE Workbench, which 
created by Dr. Aguirrcgabiria is part of the Physics Academic Software. Internally this package 
uses a Dormand-Prince method of order 5 to integrate differential equations. Further information at 
http: / / www.webassign.net /pas / ode/odewb.html 

This model has been exhaustively analyzed before, revealing a notable richness of emergent temporal 
structures which included the three main routes to chaos, as well as a multiplicity of stable coexisting 
states, see for more details [25 11321155] . 

Transfer Entropy 

TE allows for a quantification of how much the temporal evolution of the activity of one enzyme helps to 
improve the future prediction of another. The oscillatory patterns of the biochemical metabolites might 
have information which can be read-out by the TE. 

For a convenient derivation, let generally assume that each of the pairs of enzymatic activity is 
represented by the two time series X = {xt]J = i and Y = {yt}t=i ■ Here, Xt is the state value of the variable 
X in time t, and similarly for y t . Let I(X P ,Y P -> X F ) = - J2 Xt+uXt , yt p { x t+i,Xt,yi)log 2 P(x t+1 \x t} y t ) 
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be the amount of information required to predict the future of X (X F ) known both the pasts of X 



and Y (X p and Y p ). Analogously, let I(X P -> X F ) = - £ X( P(x t+1 ,x t ) log 2 P(x t+1 \x t ) be the 



amount of information required to predict the future of X known only its past. The difference I(X P — > 
X F ) — I(X P ,Y P — > X F ) is by definition the transfer entropy from Y to X, denoted by TEy_>x- It 
quantifies the amount of information in digits that Y adds to the predictability of X . 

Rewriting the conditional probabilities as the joint probability divided by its marginal, one obtains 
an explicit form for the Transfer Entropy: 



The formula (|6|) is fully equivalent to the Mutual Information between X F and Y p conditioned 
to X p . Thus, TEy^x = I(X F ,Y P \X P ), and consequently, Transfer Entropy says about how much 
information the inclusion of Y p adds to the prediction of X F only considering X p , ie. I{X F , Y P \X P ) = 
H(X F \X P ) — H(X F \X P , Y p ). Therefore, TE is fully quantifying the information flows between pairs of 
variables. The values of TE were normalized between and 1. 

It is important to remark that the TE from X to Y is different to the one from Y to X, ic. the 
effective connectivity is asymmetric, adding a directionality in time which accounts for a particular case 
of directed graphs, the graph of information flows between pairs of enzymes. 

Alternatively to the Transfer Entropy, effective connectivity can be obtained using Granger Causality 
[41) , which makes emphasis on how much from the past of one variable the predictability of its future is 
improved by adding the past of another variable. Recently, it has been proved that in the case of Gaussian 
variables both Transfer Entropy and Granger Causality are measuring exactly the same [42]. Therefore, 
the information flows based on Transfer Entropy and the Granger causality interactions coincide for 
Gaussian variables. 

Mutual Information and Redundancy 

MI quantifies how much the knowledge of one variable reduces the entropy or uncertainty of another. 
Therefore, MI says about how much information the two variables are sharing. The strongest point 
of the MI is that it extends functionality to high order statistics |28| . Its definition is MI(X,Y) = 




(6) 
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H(X) - H(X\Y), where H(X\Y) = H(X,Y) - H{Y) is the conditional entropy of X given Y. It 
accounts for the remaining uncertainty in X knowing the variable Y. We referred H(X) and H(X, Y) as 
respectively the joint and marginal (Shanon) entropies. 

For statistical independent X and Y variables one has MI(X, Y) = 0. The other limit satisfies 
MI(X, X) = H(X), because of H(X\Y) = 0. Therefore, the MI of two variables is bounded and satisfies 
that < MI(X, Y) < H(X). High values of MI mean that the redundancy in information between the 
two variables is large. The values of MI were normalized between and 1. 

Number of bins vs Statistical significance 

For all the probabilities used in both Transfer Entropy (TE) and Mutual Information (MI) we used a 
number of bins of 10. As it is well-known, the calculation of these probabilities is sensitive to the number 
of bins. Instead of tuning it as a control parameter to compute the probabilities, we preferred explored 
the statistical significance of the computed values. This was achieved by comparing both the TE and MI 
between the two series of enzymatic activity, say X and Y, with the values obtained when considering a 
random permutation of Y, what we called, the shuffled Y. The values of both TE and MI shown in Tables 
I and III were larger than those calculated in the shuffled situation (for both TE and MI, pvalue=0, for 
50 different samples). 
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Figure 1. Multi-enzyme instability-generating system of yeast glycolysis. The main 
irreversible enzymatic processes are arranged in series: Ei (hexokinase) , E2 (phosphofructokinase) and 
E3 (pyruvatekinase) . S, Pi, P[, P 2 , P' 2 and P3 denote, respectively, the concentrations of glucose, 
glucose-6-phosphatc, fructose 6-phospfate, fructose 1,6-bisphospfate, phosphoenolpyruvate and 
pyruvate, qi is the rate first-order constant for the removal of Pi; q 2 is the rate constant for the sink of 
the product P3. The model includes the feedback activation of E 2 and the feedback inhibition of E3. 
The ATP is consumed by Ei and recycled by E3. 
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Figure 2. Glycolytic route to chaos and dynamical effective connectivity. Left Panel: The 
time evolution of the E2 activity (the normalized concentration /3, fructose 1,6-bisphospfatc) shows a 
quasi-periodic route to chaos when varying the amplitude of the periodic input-flux from A= 0.001 
(top) to A= 0.023 (bottom), (a) Periodic pattern, (b) Quasi-periodic oscillations, (c) Complex 
quasi-periodic motion indicating the beginning destruction of the periodic behavior, (d) Deterministic 
chaos. All series are plotted after 10000 seconds. Right Panel: Effective connectivity of the system for 
the same values of A in the left panel. The strength of effective connectivity is plotted with arrows 
width proportional to the Transfer Entropy divided by the maximum value (red arrow), cf. results given 
in Table I. Black dashed circles at the TE from E3 and E2 emphasize that the strength of Information 
flows is not the same, but varies trough the quasi-periodic route to chaos. 
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Figure 3. Total information flows and the functional invariant. Bars represent the total 
information flow, defined per each enzyme as the total outward TE minus the total inward. For 
A=0.021 and E2 an schematic visualization of the calculation of this flow is shown (bottom graph of the 
panel). The functionality attributed for each enzyme is an invariant and preserved along the route, ie. 
E2 is a source, E3 is a sink and Ei has a quasi-zero flow. 



Tables 

Table I. Values of normalized Transfer Entropy 





From Ei 


From E 2 


From E 3 


To Ei 




[0.73;0.88;0.72;0.76] 


[0.74;0.80;0.68;0.74] 


To E 2 


[0.76;0.80;0.70;0.72] 




[0.58;0.84;0.61;0.66] 


To E 3 


[0.78; 0.86;0.74;0.75] 


[1.00;1.00;1.00;1.00] 





Table II. Values of total information flows 



Ei 


[0.06; -0.02;0.04;-0.04] 


Quasi-zero flow 


E 2 


[0.40; 0.24;0.41;0.39] 


source 


E 3 


[-0.46; -0.22;-0.45;-0.35] 


sink 
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Table III. Time Correlations 





Ei 


E2 


E 3 


Ei 


[1.00;1.00;1.00;1.00] 


[-0.65;-0.66;-0.64;-0.63] 


[-0.66;-0.66;-0.66;-0.66] 


E2 


[-0.65;-0.66;-0.64;-0.63] 


[1.00;1.00;1.00;1.00] 


[0.90;0.90;0.90;0.90] 


E3 


[-0.66;-0.66;-0.66;-0.66] 


[0.90;0.90;0.90;0.90] 


[1.00;1.00;1.00;1.00] 



Table IV. Values of normalized Mutual Information 





Ei 


E2 


E3 


Ei 


[1.00;1.00;1.00;1.00] 


[0.52;0.49;0.45;0.44] 


[0.48;0.49;0.45;0.44] 


E2 


[0.52;0.49;0.45;0.44] 


[0.85;0.84;0.85;0.86] 


[0.47;0.46;0.45;0.45] 


E3 


[0.48;0.49;0.45;0.44] 


[0.47;0.46;0.45;0.45] 


[0.76;0.74;0.76;0.78] 



